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Abstract 

We study the production of a prompt J/ijj meson in association with a prompt 
photon in ep deep-inelastic scattering within the factorisation formalism of non- 
relativistic quantum chromo dynamics (NRQCD) and demonstrate that this process 
provides a clean probe of the colour-octet mechanism at DESY HERA. Our analysis 
is based on an updated set of non-perturbative NRQCD matrix elements obtained 
through a joint fit to data on charmonium inclusive hadroproduction from runs I 
and II at the Fermilab Tevatron. 
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1 Introduction 



Since the discovery of the J/if> meson in 1974, charmonium has provided a useful labora- 
tory for quantitative tests of quantum chromodynamics (QCD) and, in particular, of the 
interplay of perturbative and non-perturbative phenomena. The factorisation formalism 

of nonrelativistic QCD (NRQCD) [2] provides a rigorous theoretical framework for the 
description of heavy- quarkonium production and decay. This formalism implies a separa- 
tion of short-distance coefficients, which can be calculated perturbatively as expansions in 
the strong-coupling constant a s , from long-distance matrix elements (MEs), which must 
be extracted from experiment. The relative importance of the latter can be estimated by 
means of velocity scaling rules; i.e., the MEs are predicted to scale with a definite power of 
the heavy-quark (Q) velocity v in the limit d<1. In this way, the theoretical predictions 
are organised as double expansions in a s and v. A crucial feature of this formalism is that 
it takes into account the complete structure of the QQ Fock space, which is spanned by 
the states n = 2S + 1 l / y> with definite spin S, orbital angular momentum L, total angular 
momentum J, and colour multiplicity a = 1,8. The hierarchy of the MEs predicted by 
the velocity scaling rules is explained for the J/if>, ip', and Xcj mesons in Tabled In 
particular, this formalism predicts the existence of colour-octet (CO) processes in nature. 
This means that QQ pairs are produced at short distances in CO states and subsequently 
evolve into physical, colour-singlet (CS) quarkonia by the non-perturbative emission of 
soft gluons. In the limit v — *■ 0, the traditional CS model (CSM) [3J is recovered in 
the case of .S-wave quarkonia. The greatest triumph of this formalism was that it was 
able to correctly describe |4l5j the cross section of inclusive charmonium hadroproduction 
measured in pp collisions at the Fermilab Tevatron ;6 3 , which had turned out to be more 
than one order of magnitude in excess of the theoretical prediction based on the CSM. 
Apart from this phenomenological drawback, the CSM also suffers from severe concep- 
tual problems indicating that it is incomplete. These include the presence of logarithmic 
infrared singularities in the 0(a s ) corrections to P-wave decays to light hadrons and in 
the relativistic corrections to S"-wave annihilation [7] , and the lack of a general argument 
for its validity in higher orders of perturbation theory. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and vN deep-inelastic 
scattering (DIS), e + e~ annihilation in the continuum, Z-boson decays, 77 collisions, and 
6-hadron decays may be found in the literature; for reviews, see jHJ. Furthermore, the 
polarisation of if}' mesons produced directly and of J/ip mesons produced promptly, i.e., 
either directly or via the feed-down from heavier charmonia, which also provides a sensi- 
tive probe of CO processes, was investigated |9|lUj . Until recently, none of these studies 
was able to prove or disprove the NRQCD factorisation hypothesis. However, HI data of 
e + p — > e + J/ip + X in DIS at the DESY Hadron Electron Ring Accelerator (HERA) 
[TTrb] and DELPHI data of 7 + 7 -> J/ip + X at the CERN Large Electron Positron Col- 
lider (LEP2) provide first independent evidence for it by agreeing with the respective 
NRQCD predictions [TliTH| . 



2 



Table 1: Values of k in (<D H [n]\ oc v k for H = J ftp, if)', x. 
k J/ip, ip' Xcj 
3 3 S( 1} 

5 3 Pj 1} , 3 ^ 8) 

7 lo(8) 3c(8) 3p(8) 



In this paper, we identify the DIS process 

e + p^ J/ip + -? + X (i; 

as a clean probe of the CO mechanism and propose its experimental study at HERA II 
In fact, among the partonic subprocesses contributing at LO, 

W + 7, 



e + 7 
e + g 



e + cc 
e + cc 



(2) 
(3) 



the latter is by far dominant because, in the relevant range of x and Q 2 , the density of 
gluons in the proton is so much higher than the one of photons that the 0(v 4 ) suppression 
of [o J l^ 3 s{ 8 ^ ^ relative to (o J ^ 3 S[^ ^ (see Table [TJ is inconsequential. The emission 
of photons off the proton can happen either elastically or inelastically, i.e., the proton 
stays intact or is destroyed. In both cases, the parton density function (PDF) can be 
evaluated in the Weizsacker- Williams approximation |16|17j . Besides electromagnetic 
proton interaction, also diffract ive scattering off the proton, via pomeron exchange, allows 
for CS processes. However, such events will be accumulated at the border of the phase 
space, at z<^l, where z is the inelasticity variable defined in Sect. By the same 
token, they can be eliminated from the experimental data set by applying an appropriate 
acceptance cut on z. In this paper, we assume this to be done. 

The potential of J/ifj plus photon associated production to probe the CO mechanism 
was already investigated for photoproduction in ep scattering jTH] . In that case, however, 
the bulk of the cross section, more than 2/3 for J /if) transverse momenta pt > 1 GeV (see 



cc 



3 S[^ +7 in resolved photoproduction. 



cc 



3 S[ 8 ^ +7, only makes 



Table II of [HI), is due to the CS channel g+g 
Specifically, the photoproduction analogue of process (jSJ), 1+g 
up 1/5 of the total cross section ^H]. On the other hand, the probability of a photon to 
appear resolved rapidly decreases with its size l/Q 2 , so that the situation encountered in 
[TH] is subject to a dramatic change as one passes from photoproduction to DIS. 

In PH], prompt J/ if) plus photon associated photoproduction in 77 collisions was 
studied in next-to-leading order (NLO), and sizeable corrections to the cross section were 
found. Since the present analysis is at an exploratory level, considering process (JTJ for 
the first time, we stay at LO. A NLO analysis of process (JTJ) would be more involved than 
in the case of jTH], due to the presence of Q 2 as an additional mass scale, and is left for 
future work. 
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The paper is organised as follows. In Sect. |2l we collect the formulas from which the 
LO cross section of process (JTJ) can be evaluated. In Sect. El we update the extraction 
of the NRQCD MEs of the J/ip, ip', and Xcj mesons in j2j by including in the fit CDF 
data from Tevatron run II (2Hj besides that from run I In Sect. HI we then present our 
predictions for the cross section of process (pQ) under HERA II experimental conditions 
and demonstrate that this is an excellent probe of the CO mechanism. Our conclusions 
are contained in Sect. El 



2 Analytic results 

We now present our analytic results for the cross section of process $IJ). We work at 
LO in the parton model of QCD with n/ = 3 active quark flavours and employ the 
NRQCD factorisation formalism pQ to describe the formation of the J/ip meson. We 
start by defining the kinematics. We denote the four-momenta of the incoming lepton 
and proton and the outgoing lepton, J/ijj meson, and photon by k, P, k', p^, and p', 
respectively. The parton a struck by the virtual photon (7*) carries four-momentum 
p = xP. The virtual photon has four-momentum q = k — k', and it is customary to 
define Q 2 = —q 2 > 0, y = q ■ P/k ■ P, and z — p^ ■ P/q • P. In the proton rest 
frame, y and z measure the relative lepton energy loss and the fraction of the virtual- 
photon energy transferred to the J/ip meson, respectively. We neglect the masses of 
the proton, lepton, and light quarks, call the one of the J/ip meson M^,, and take the 
charm-quark mass to be m c = In our approximation, the proton remnant X has 

zero invariant mass, M\ = [P — p) 2 = 0. The centre-of-mass (cm.) energy squares of 
the ep and 7*p collision are S = [k + P) 2 and W 2 — (q + P) 2 = yS — Q 2 , respectively. 
As usual, we define the partonic Mandelstam variables as s = (q + p) 2 = xyS — Q 2 , 
t — (q — pip) 2 = —xy{l — z)S, and u = (p — p^) 2 = Ml — xyzS. By four-momentum 
conservation, we have s + t + u = M 2 — Q 2 . In the 7*p CM frame, the J/ip meson has 
transverse momentum and rapidity 

^t(su + Q 2 M 2 ) 

Vt = —^q 2 — ' (4) 

y* h = -ln—, ^—^ '- + -ln , (5) 

* 2 s(M 2 -t)+Q 2 M 2 2 s y J 

respectively. Here and in the following, we denote the quantities referring to the 7*p CM 
frame by an asterisk. The second term on the right-hand side of © originates from the 
Lorentz boost from the 7*0 CM frame to the 7*p one. Here, y^ is taken to be positive 
in the direction of the three-momentum of the virtual photon, in accordance with HERA 
conventions |ll|12j . 

The cross sections of processes (j2J) and (jHJ) may be conveniently calculated by applying 
the covariant-projector method of They are related to the one of e+g — > e+cc 3 S[ 8 ^ + 
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g, given in (13) of [Hj, by 



d 3 a 



(e + 7 



dy dQ 2 dt 

64 a 2 d 3 a 



e + cc 



9 a 2 dy dQ 2 dt 

d 3 a ( 

[e + g-. 



(e + g — > e + cc 



+ 3 



<iQ 2 (it 
„ a d 3 cr 



a s dy dQ 2 dt 



-> e + cc 



(6) 



where the proportionality factors account for colour and coupling adjustments. 

According to the factorisation theorems of the parton model and NRQCD, the cross 
section of process ((TJ) is then evaluated as 



d 2 



a 



dy dQ 2 
fi 



[e + p -> e + J/V> + 7 + X) 

r-0 



dx 



(Q 2 +A€)/(2/S) 



dt 



(s+Qi)(s-Ml)/s 



E/-/p(^^)E(° j/v W 



d 3 



x 



dQ 2 dt 



(e + a — > e + cc[n] + 7). 



(7) 



where the sums run over (a,n) = (7, 3 Si ) , (g, 3 Si\ f a / p (x,M) is the PDF of parton 
a in the proton at factorisation scale M, and (d 3 a/dy dQ 2 dt) (e + a — > e + cc[n] + 7) 
are given by (JHJ). The kinematically allowed ranges of y and Q 2 are M 2 /S < y < 1 and 
< Q 2 < - M|, respectively 

Prompt J/^ production may be conveniently described by inserting in (J7J) the effective 
MEs specified in (20) of [Hj. 



3 Determination of the MEs 

The recent CDF measurement of prompt J /if) inclusive hadroproduction in run II at the 
Tevatron (with \/S = 1.96 TeV and \y^\ < 0.6) [2U] allows us to update and improve 
our knowledge of the CO MEs. Previous determinations by one of us in collaboration 
with Kramer [22] and with Braaten and Lee [5] were only based on run I data (with 
y/S = 1.8 TeV and \y^\ < 0.6) 0. However, the latter data were more detailed because 
the prompt J /if) sample was explicitly broken down into direct J /if) mesons, feed-down 
from if)' mesons, and feed-down from XcJ mesons. In order to make maximum use of the 
available information, we perform a joint fit to the data from runs I and II, which otherwise 
proceeds along the lines of 0. First, the CS MEs (0'^ nS) \ 3 S{ 1 A), with n = 1,2, and 
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and Xc2 

Qip(nS) 



3 P f 



(1) 



are extracted from the measured partial decay widths of ip(nS) — > / + + I 



7 + 7 



3p(8) 
-in 



and (0* 



c0 [3g-( 



(8) 



respectively. Then, the CO MEs (0^ nS ^ 3 S{ 8) ), (0^ nS ^ 



15.(8) 



are fitted to the p? distributions of ip(nS) and Xcj 
inclusive hadroproduction |6I20| and the cross-section ratio & Xc2 /a Xcl |2I] measured at the 
Tevatron. In contrast to the run I data jHj, the run II data [20] reach down to px = 0. 
It turns out that the run II data in the newly covered low-p^ range are comparable to or 
below the CSM prediction, so that their inclusion would spoil the fit. Therefore, we on 
include the 14 data points with > 4.25 GeV. Since the fit results for (p^ nS) ' 

and /(9^( nS ) 3 Pq ) are strongly correlated, we consider the linear combination 



15.(8) 



j^ip(nS) _ /Qip(nS) 



15,(8) 



ml \ 



3 P, 



(8) 



where the value of r is chosen so that the error on Mf^ nS ^ is minimised. The minimisation 
in x 2 is performed exactly, by solving a set of seven linear equations for the seven unknown 
CO MEs. 

The relevant partonic cross sections may be found in j^j. We take the charm-quark 
mass to be m c = (1.5 ± 0.1) GeV and adopt the relevant feed-down and leptonic decay 
branching fractions from j2H|. As for the proton PDFs, we use the latest LO sets by Mar- 
tin, Roberts, Stirling, and Thorne (MRST2001LO) and the Coordinated Theoretical- 
Experimental Project on QCD 
(CTEQ6L1) [26 a . For consistency, we employ the one-loop formula for a^(/j,) and choose 
the asymptotic scale parameter to be Aq^ d = 253 MeV (246 MeV), appropriate for the 
MRST2001LO (CTEQ6L1) PDFs. We identify the renormalisation scale /x and the fac- 
torisation scale M with the charmonium transverse mass my = \jp\ + M^. 

The sets of MEs thus obtained with the MRST2001LO and CTEQ6L1 PDFs are 
summarised in Table El The respective values of x 2 P er degree of freedom are 49/52 
and 50/52. The quoted errors are of experimental origin only. Comparing Tabled with 
Table I of [Oj, where previous MRST and CTEQ PDFs were used, we observe that the 



MEs obtained with CTEQ PD 



MRST values of 
than in jU]. 



Qip(nS) 35.(8) 



's are only moderately changed. On the other hand, the 
and (O Xc0 3 S[ 8 ^ \ in Table 121 are appreciably smaller 



4 Prompt J /if) plus photon associated electroproduc- 
tion 

We are now in a position to present our theoretical predictions for the cross section of 
process (0) under HERA II experimental conditions. They are evaluated from the formulas 
listed in Sect. El with the inputs and conventions described in Sect. except that we now 
set p, — M — ^\fQ 2 + M 2 and vary the scale parameter £ between 1/2 and 2 about the 
default value 1. Since (JJJ) is sensitive to a different linear combination of (o^^ nS ^ 
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Table 2: NRQCD MEs of the J/ip, ip', and XcJ obtained as described in the text using 
the MRST2001LQ [25] and CTEQ6L1 gH] PDFs. The errors are exp erimental only 

CTEQ6L1 




MRST2001LO 



1.5 ±0.1 GeV 3 



1.4 ±0.1 GeV 3 



(3.9 ± 1.4) 


X 


10" 


- 4 GeV 3 


(2.3 ±0.2) 


X 


io- 


- 3 GeV 3 


(6.0 ±0.1) 


X 


io- 


2 GeV 3 


(7.3 ±0.2) 


X 


io- 


- 2 GeV 3 


(6.7 ±0.5) 


X 


io- 


- 1 GeV 3 


(6.7 ±0.5) 


X 


io- 


- 1 GeV 3 


(8.8 ± 1.2) 


X 


io- 


4 GeV 3 


(2.0 ±0.2) 


X 


io- 


- 3 GeV 3 


(1.1 ±0.1) 


X 


io- 


2 GeV 3 


(1.0 ±0.2) 


X 


io- 


2 GeV 3 


(1.2 ±0.1) 


X 


10" 


- 1 GeV 5 


(1.2 ±0.1) 


X 


io- 


- 1 GeV 5 


(4.7 ±0.6) 


X 


10" 


4 GeV 3 


(1.1 ±0.1) 


X 


io- 


- 3 GeV 3 



and (0^ nS ^ ^Pq ) than the one appearing in (jHJ), we write 



Qip(nS) 
Qip(nS) 



"lo(8) 
3p(8) 



KMf {nS \ 



r 



(9) 



and vary k between and 1 around the default value 1/2. 

In order to estimate the theoretical uncertainties in our predictions, we vary the un- 
physical parameters £ and k as indicated above, take into account the experimental errors 
on m c and the default MEs, and switch from the MRST2001LO PDFs |25j . which we take 
as our default, to the CTEQ6L1 ones [2E], properly adjusting Aq^d and the MEs. We 
then combine the individual shifts in quadrature, allowing for the upper and lower half 
errors to be different. 

In the laboratory frame of HERA II, electrons or positrons with energy E e = 27.5 GeV 
collide with protons of energy E p 



920 GeV, yielding a cm. energy of V S = 2 s jE e E p = 
318 GeV. For defmiteness, we adopt the experimental acceptance cuts from the recent 
HI analysis of inclusive J/ij) electroproduction [TT], which include 2 < Q 2 < 100 GeV, 
50 < W < 225 GeV, 0.3 < z < 0.9, and pf > 1 GeV 2 . We consider cross section 
distributions in Q 2 , W, z, p^ 2 , yl, pj,, and y^, where the last four variables refer to the 
J/ip meson, adopting the binning from [TT] . 

Our results are displayed in Fig. ^ where the NRQCD predictions (solid lines) are 
compared with the CSM ones (dashed lines). In each case, the theoretical uncertainties are 
indicated by the shaded bands. As expected from our discussion in Sect. UJ the NRQCD 
predictions vastly exceed the CSM ones, by almost one order of magnitude, The gaps 
are considerably larger than the theoretical errors. On the other hand, the shapes of the 
various cross section distributions come out very similar in both approaches, which may 
be understood by observing that the partonic cross sections of the CS and CO channels in 
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(J!))) are proportional to each other. From Fig. ^ we read off that the integrated NRQCD 
cross section is of order 10 fb. Given that the integrated luminosity to be collected by 
the end of HERA operation amounts to about 1 fb -1 , one thus expects about 10 signal 
events. 

5 Conclusions 

We studied the electroproduction of prompt J ftp mesons in association with prompt 
photons in ep collisions under HERA II kinematic conditions to LO in the NRQCD [2\ 
factorisation formalism 1\. We considered cross section distributions in all variables 
of current interest |llll2j . including Q 2 , W, z, pj?, y^, p\, and y^. As input for our 
calculation, we used updated information on the NRQCD MEs extracted from a combined 
fit to data on inclusive charmonium hadroproduction collected by the CDF Collaboration 
in runs I [H] and II (2H| at the Tevatron. 

As a result of our study, we could identify prompt J/ip plus photon associated elec- 
troproduction as a useful probe of the CO mechanism. In fact the NRQCD predictions 
turned out to exceed the CSM ones by almost one order of magnitude. Unfortunately, the 
cross section distributions in both theories have very similar shapes, so that the NRQCD 
to CSM ratio cannot be further enhanced by specific acceptance cuts. Should this pro- 
duction process be experimentally observed with the rate predicted by NRQCD, then 
this would provide strong evidence in favour of the existence of CO processes in nature. 
In view of the moderate integrated cross section of order 10 fb, this is a challenging en- 
deavour. This is an example of a study that would benefit from the extension of HERA 
operation beyond the summer of 2007. 
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